Discrete surface solitons in two-dimensional anisotropic photonic lattices 
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We study nonlinear surface modes in two-dimensional anisotropic periodic photonic lattices and 
demonstrate that, in a sharp contrast to one-dimensional discrete surface solitons, the mode thresh- 
old power is lower at the surface, and two-dimensional discrete solitons can be generated easier near 
the lattice corners and edges. We analyze the crossover between effectively one- and two-dimensional 
regimes of the surface-mediated beam localization in the lattice. 

PACS numbers: 



Surface modes have been studied in different branches 
of physics; in guided wave optics surface states were pre- 
dicted to exist at interfaces separating periodic and ho- 
mogeneous dielectric media 0. The interest in studying 
surface waves has been renewed recently because the in- 
terplay of discreteness and nonlinearity can facilitate the 
formation of discrete surface solitons 0, at the edge of 
the waveguide array. That can be understood as the lo- 
calization of a discrete optical soliton near the surface Q 
for powers exceeding a certain threshold value, for which 
the repulsive effect of the surface is balanced. A similar 
effect of light localization near the edge of the waveg- 
uide array and the formation of surface gap solitons have 
beenpredicted and observed for defocusing nonlinear me- 
dia [IQ. 

It is important to analyze how the properties of nonlin- 
ear surface waves are modified by the lattice dimension- 
ality, and the first studies of different types of discrete 
surface solitons in two-dimensional lattices 0, El E3 
revealed, in particular, that the presence of a surface in- 
creases the stability region for two-dimensional (2D) dis- 
crete solitons ^3 an( i t ne threshold power for the edge 
surface state is slightly higher than that for the corner 
soliton 0- 

In this Letter we consider anisotropic semi-infinite two- 
dimensional photonic lattices and study the crossover be- 
tween one- and two-dimensional surface solitons empha- 
sizing the crucial effect of the lattice dimensionality on 
the formation of surface solitons. 

We consider a semi-infinite 2D lattice [shown schemat- 
ically in FigEIa) below], described by the system of 
coupled-mode equations for the normalized amplitudes 

Un,m PII2, 




FIG. 1: Examples of (a) corner, (b,c) two edge surface, and 
(d) central localized nonlinear modes for A = 5. 



fine the lattice coupling as follows: 



u 2 , m n = 1, m > 1, 

^n+l,m d~ V"n,—l 7 m Ti ^> \, 

U n .2 Tl > 1, m = 1, 

(u n ,m+l + U n ,m-l) Wl > I, 



i d""' T " + (y„ + V m ) u n m + \u n , m \ 2 u n , m = 0, (1) 
where £ is the normalized propagation distance. We de- 



where a characterizes the lattice anisotropy. 

Linear lattice waves of the form u„. m (£) = 
Mo sin(fcri) sin(qm) exp(i/3£) satisfy the dispersion rela- 
tion /3kq = 2(cosfc + acosq). In the nonlinear case, 
we look for localized stationary solutions of the form 
(0 — «n,m cxp(iA£), where the amplitudes 

^n,m are 

real, and A is the nonlinear propagation constant. For a 
given A, localized solutions are found in a 15 x 15 lattice 
by using the Newton- Raphson method. 

We calculate the power threshold Pth that character- 
izes the discrete solitons in 2D lattices [3. We study 
three different modes: corner [Fig. IHa)], edge surface 
[Fig. mb,c)] and central [Fig.^d)] localized modes. The 
corner and edge modes represent 2D surface localized 
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modes, and the central mode corresponds to a 2D dis- 
crete soliton of an infinite lattice. We find the total power 
P = J2 n m \ u n,m\ 2 of all those modes and perform a lin- 
ear stability analysis 0, 0| for each solution [Fig. Ufa)] . 

We observe that the threshold power P t h for the sur- 
face modes is smaller than the power corresponding to 
the central mode, with the corner state having the small- 
est P t h- This interesting feature was previously ob- 
served for a single nonlinear impurity placed near a 
boundary of a 2D lattice 0], and it appears also in the 
anisotropic model. Therefore, in a sharp contrast with 
one-dimensional (ID) surface solitons, the surface of a 
2D lattice creates an effectively attractive potential for 
the localized modes that reduces the threshold power for 
the mode localization. 




FIG. 2: (Color online) (a) Power diagram of the central (con- 
tinuous), edges surface (dashed), and corner (dotted) local- 
ized modes for the isotropic case (a = 1). Inset: Schematic of 
a semi-infinite two-dimensional lattice. (b,c) Examples of the 
threshold power landscapes P t h(n,m) for a = 1 and a = 0, 
respectively. 

The linear stability analysis of the surface localized 
modes coincides with the Vakhitov-Kolokolov stability 
criterion [l^: the surface modes are stable for dP/dX > 
0, and unstable otherwise. To study in more details how 
the threshold power P t h varies for the modes localized 
at different points of the lattice, we compute the power 
threshold landscape: We look numerically for one-peak 
localized modes centered at different sites of one-quarter 
of the whole lattice and find the first stable solution at 
the threshold power which corresponds to Pth at the site 
(n, m). In Fig.[2JF>,c), we show the examples of two land- 
scapes for two different values of the lattice anisotropy 
a. For a = 1, Pth grows from the minimum P t h ~ 4.9 



at the corner (1,1) to two maxima corresponding to the 
edge surface modes (8, 1) and (1, 8) and the central mode 
(8, 8). For a = 0, we observe that the required power to 
excite a localized mode decreases from the surface (n = 1) 
to the center (n = 8). The value P t /,(l,m) as 3.2 corre- 
sponds to the threshold power of discrete surface solitons 
in a semi-infinite array J4j. For ID localized modes there 
exists no power threshold in the continuum limit [l3l | , but 
in our system P t h(8,m) ~ 0.4 due to finite-size effects. 

Next, we study the effect of the lattice anisotropy on 
the power threshold of discrete surface solitons. In Fig. [21 
we show P t h for different values of the parameter a (0 < 
ex < 1) , for the four types of localized modes shown in 
Figs.DIa-d). 
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FIG. 3: (Color online) Threshold power P t h vs. a for the 
central (continuous thick), two edge surface (continuous thin 
and dotted), and corner (dashed line) localized modes. Inset: 
Pth(n,m) for a = 0.3. 

When the anisotropy parameter a grows, the threshold 
power grows for all localized modes. It means that, by 
increasing the effective dimensionality of the system, the 
required power to excite a nonlinear discrete localized 
state also increases. The threshold power for two edge 
surface modes decreases as a decreases to a « 0.3, when 
the power curves cross and diverge. For a = 0, the edge 
soliton centered at (8, 1) [Fig. QJc), dotted line in Fig. [3] 
only interacts with an effective ID array (m = 1). On 
the contrary, P t h for the edge surface mode centered at 
(1, 8) [Fig.njb), continuous-thin line in Fig.|3] approaches 
the same limit as the corner mode. For a = 0, this mode 
only interacts with an effective ID chain (m = 8) and the 
surface. An interesting feature of Fig. |3|is that it shows 
the existence of a critical value of the lattice anisotropy, 
a « 0.3 where all powers thresholds almost coincide. In 
other words, there exists a critical value of the lattice 
anisotropy where we observe a crossover between ID and 
2D lattices. 

To study this interesting effect in more details, we cal- 
culate the power threshold landscape P t h for a = 0.3 [see 
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inset in Fig. [3] , and observe that the localized mode cen- 
tered at (2, to) possesses the lowest value of P t h in the 
array and, in particular, the localized mode centered at 
(2, 1) is the mode with the lowest threshold power. For 
a < 0.3, the minimum P t h will be located outside of the 
array boundaries, and intermediate localized states (lo- 
cated between the center and the surfaces of the array) 
will be the easier states to excite. For ID semi-infinite 
arrays Q , it was shown that the power threshold for sur- 
face modes decreases as the mode moves away from the 
surface [see also Fig.^c)]. 



(a) ? = f = 3 




(b) f = f = 500 




FIG. 4: Generation of surface modes by one-site excitation 
(a = 1) at (a) center and (b) corner of the lattice. 

Finally, we study the dynamic generation of localized 



modes in two-dimensional lattices. We simulate numeri- 
cally Eqs. by using a one-site excitation u o~ n ,ri&m,m' ■ 
In Fig. 0J we show our results by plotting the evolution of 
the mode amplitude |u n , m | for two different values of £. 
As an example, we show the excitation of the corner and 
central localized modes for the initial amplitude uq = 2.5. 
If we try to excite a localized mode at the middle of the 
lattice for this value of uo, this localized state decays 
quickly, as shown in Fig. Ufa), an d the beam diffracts 
because the input power is below the effective dynami- 
cal power threshold. On the contrary, by using the same 
input power we are able to excite the corner soliton [see 
Fig. Etb)] which is dynamically stable in the propaga- 
tion. It is worth mentioning that, in the one-dimensional 
limit (a = 0) the situation becomes completely differ- 
ent: If we try to excite a surface mode at n = 1 with 
the power smaller than a threshold value, but still large 
enough, the input beam will not decay but will instead 
move to the neighboring lattice site creating a localized 
mode there . 

In conclusion, we have analyzed the properties of dis- 
crete surface solitons in two-dimensional anisotropic pho- 
tonic lattices, and studied the crossover between the one- 
and two-dimensional regimes of the surface-mediated 
beam localization in the lattice. In particular, unlike 
one-dimensional discrete surface solitons, the threshold 
power of the two-dimensional discrete soliton is lowered 
by the surface, so that two-dimensional solitons can be 
generated easier near the lattice corners and edges. 

This work was supported by Fondecyt grants 1050193 
and 7050173 and the Australian Research Council. 
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